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G\ ' Abstract 

We examine theoretically the effects of random topographical substrates on the motion of two- 

'Ch' 

dimensional droplets via appropriate statistical approaches. Different random substrate families are 

i 

and in the long-time limit is deduced and the droplet footprint is found to be a normal random 
variable at all times. It is shown that substrate roughness decreases droplet wetting, illustrating 
£>y also the tendency of the droplet to slide without spreading as equilibrium is approached. Our 

1 Ph ' theoretical predictions are verified by numerical experiments. 
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Front propagation in heterogeneous media occurs in a wide variety of areas in physics, 
ranging from transport phenomena in porous media and reaction-diffusion-advection sys- 
tems to crack propagation due to lattice defects Heterogeneities are always present in 
natural environments but quite often their precise form is unknown. In such cases they can 
be modeled as random noise signals, a more physical and practical assumption than e.g. 
periodic ones. An example of front propagation central to interfacial hydrodynamics is that 
of a moving contact line during liquid spreading on a solid substrate, where heterogeneities 
usually originate from substrate defects, either chemical {2} or topographical 

It is a fundamental problem to understand how random heterogeneities influence the 
spreading dynamics and the characteristics of contact line propagation, e.g. speed and 
location. Experimental studies on droplet spreading - a simple prototype for the study of 
contact line motion - suggest that substrates having highly irregular micro-scale features, 
commonly called rough, can influence the dynamics significantly (4]. Several theoretical 
works focused on equilibrium configurations and deterministic substrates: the study by 
Wenzel j^j on droplet equilibria on rough substrates obtained an effective contact angle 
based on simplistic energy arguments - a formula which was supported by experiments in a 
regime where the contact angles are not small j^]; in the contact line of a fluid wedge in the 
presence of a single, localized defect was determined- in [?| droplet equilibria on substrates 
with regular, periodic features were considered; Cox 8| investigated multiple equilibria of an 
infinite fluid wedge on a general non-periodic rough surface. Introducing randomness in the 
variations of the substrates is clearly a realistic way to represent roughness. There are a few 
studies in this direction but are based on phenomenological modeling ideas as for example 9] 
which regarded rough substrate patches as contact angle point sources or [lO] which looked 
into the effects of surface roughness using an Ising model. Studies on contact line dynamics 
on random substrates were also based on phenomenological ideas, e.g. where postulated 
equations describing the dynamics were utilized. Hence, to date a systematic fluid dynamics 
treatment based on rational statistical approaches is still lacking. As a result, contact line 
motion on heterogeneous substrates is far from being well understood, unlike other problems 
in continuum mechanics, such as porous media, which not only have been analyzed from first 
principles but have also motivated new mathematical approaches (cf. homogenization 12]). 

In this Letter we report the first detailed and systematic study of the qualitative effects 
of random, small-scale spatial heterogeneities on droplet motion, through the development 
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of appropriate statistical methodologies. The starting point is the recent work in [13j on the 
motion of two-dimensional (2D), partially wetting droplets over deterministic, spatially het- 
erogeneous substrates. The restriction to 2D implies that there are no transverse variations 
and the contact line is essentially treated as a set of two points. This simplified problem, 
albeit more difficult to study experimentally, is a first step towards understanding the in- 
fluence of spatial heterogeneities, but the required analysis still remains highly nontrivial. 
The model for the droplet motion was obtained rigorously from the governing hydrodynamic 
equations. Assuming small contact angles, a long-wave expansion in the Stokes flow regime, 
yields a single equation for the evolution of the droplet thickness H (x, t) over a substrate 
i] (x) . In dimensionless form the equation reads: 

d t H + d x [H 2 (H + A) dl (H + v )] = 0, (1) 

where A <C 1 is the non-dimensional slip length originating from the Navier model imposed 



to alleviate the stress singularity that occurs at the moving contact line [I4j . Equation (JTJ) 
describes the capillarity driven spreading that is resisted by fluid viscosity: the term d t H 
results from the viscous fluid motion and the term d x [-} represents the effect of surface tension 
and also accounts for the substrate curvature which contributes an additional capillary 
pressure, —d^rj. The spatial coordinate x, and time t, are made non-dimensional by the 



characteristic length L = a/ ~Aj (2a s ) and time r = 3/j,L/( r ya^), respectively, where A is the 
droplet cross-sectional area, a s is the equilibrium angle prescribed by Young's law, \i is the 
viscosity and 7 is the surface tension. r\{x) and H(x) are scaled by La s and A by La s /3 
respectively. Equation (JTJ is subject to a constant volume constraint and the boundary 
conditions at the contact points that the droplet thickness vanishes and the angle the free 
surface makes with the substrate remains equal to its static value, a s . For quasi-static 
spreading, a singular perturbation method was employed in [3] to asymptotically match 
the solution in the bulk of the fluid with the solution in the vicinity of the contact lines. 
This led to a set of two integro-differential equations (IDEs) for the time-evolution of the 
right and left contact points at x = a± (t), respectively. Their validity has been confirmed 



by comparisons with numerical solutions to the full equation in (CQ) 131 ]. 

Here we use the set of IDEs to investigate the case where rj (x) is a random function. We 
take 1 77 (x) I Cl assuming also that its variations occur at lengthscales that are much longer 
than A. The requirement of substrate smoothness together with the fact that the 'noise' 
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FIG. 1: (a) Sample substrate realizations using ([2]) for ?7o = 10" 3 and k = 10, 20 and 40. (b) 
Top plot. Solid line: experimental profile from Hitchcock el al. |3[] for an alumina sample. Dashed 
line: approximation obtained by projecting the experimental profile onto N = 26 harmonics in ([2]). 
Lower plot: sample substrate realization for the same rjo and ko, but N = 1000. 

is spatial and enters the equations in a nonlinear fashion, precludes casting the problem 
into the standard Langevin-Fokker- Planck formalism, which is usually employed to study 
randomly perturbed dynamical systems flo| . 

The primary fundamental difficulty with the substrate is the development of a random 
representation that can have a large frequency content and at the same time it is differen- 
tiable. A convenient representation is the random function, 

= ~/W ^ am sm ~W X cos ~W X ' 

* m=l 

where ^ 1 an d fc are the characteristic amplitude and wavenumber of the substrate, re- 
spectively, and iV is a large positive integer. Here the coefficients a m and f3 m are statistically 
independent normal variables with unit variance. It is readily seen that rj (x) is a periodic 
function with period 27rN/ko, but we eventually take iV — > oo so that this periodicity is 
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lost and r](x) approaches a band-limited white noise (see Fig. [Tk; in this limit, continuity 
of (j2J) and all its derivatives can be carefully shown by Kolmogorov's continuity theorem). 
Representations of this form, are also invoked to look into noise effects in other contexts as 
they occur, for example, in electrical current signals or in black-body radiation [l6|. 

An attractive feature of the stochastic representation in (j2J) is that it generates an infinite 
family of substrate realizations parameterized by two parameters, 770 and ko, which are 
often reported in experimental studies when characterizing a rough substrate (in Fig. [lb, 
r]o ~ 0.78 — 2.1/zm, wavelength l ~ 14 — 42//m). These two parameters may be also be used 
to compute the 'roughness coefficient of the substrate', r, defined as the mean of the ratio of 
the actual surface area over its projected area. In the limit of f]oko <C 1 we find r = l+ri^kn/Q, 
to be contrasted with r = I+^q/Cq/ '4 obtained with the pure harmonic 77 (x) = r/ cosk x |l7 |. 
Comparison with experimental substrate profiles determined by Hitchcock et al. j^] shows 
that (j5J) can be used to realistically represent actual rough substrates. For example, given the 
experimental profile in Fig. [Ib (solid line of the upper plots), ko is readily determined from 
ko = 2ixn/ a/5/3 and 770 from 770 = \/Jjf), where n corresponds to the number of maxima 



per unit length [16!] and (•) denotes an ensemble average over all substrate realizations. To 
obtain the dashed profile as an approximation to the experimental profile using (jSJ), the finite 
length of the profile is matched to the period of @ to get iV = 26, which then allows us to 
determine the a m and f3 m by projecting the experimental profile onto their corresponding 
harmonics. On the other hand, the lower plot in Fig. [lb is generated with random a m and 
f3 m with N = 1000, where 770 and ko are kept the same so that the lower plot belongs to the 
same substrate family with the upper one. For the statistical analysis that follows, a large 
number of substrate realizations, typically 20, 000, will be utilized. 

Different substrate descriptions might have been used that have been observed in pro- 
filometry measurements, as for example representations that exhibit statistical self-affinity, 
whose spectral density follows the power law oc k 2D ~ 5 [2^], where 1 < D < 2 is the frac- 
tal dimension. The differentiability requirement of the substrate representation together 
with the fact that in reality a self-afnne structure cannot persist for all lengthscales requires 
imposing lower and upper wavenumber cutoffs. However, as we shall see later the results 
for the statistics of the contact lime motion are qualitatively the same regardless of the 
representation we employ. 

To facilitate the analysis, we introduce the contact line fluctuation e and droplet shift I 
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FIG. 2: Droplet lying between a_ < x < a+ on a random topographical substrate. The droplet 
shift, £ = ^(a+ + a_), is the distance the droplet midpoint moves away from x = and the contact 
line fluctuation, e = ^(a+ — ci-) — xq, measures deviations of the contact line location away from 
the flat-substrate radius, xo(t). 

along the substrate (see Fig. [2]), defined in terms of the contact line locations as 

e = \ ( a + ~~ a ) ~ x o an d t — | («+ + 0-) , (3) 

where xq is the contact line location when the droplet spreads on a flat substrate symmet- 
rically about x = 0, which approaches v^3 in the long time limit. The droplet shift is a 
measure of the distance the droplet midpoint is displaced from x = 0, whereas the contact 
line fluctuation measures the deviation of the droplet radius from the flat-substrate radius. 
As the amplitude of the topographical features is taken to be small, we also expect e <C 1. 
By linearizing about the flat-substrate equilibrium we obtain 

e = —jL= y~] ( a ™ sin X m £ + (3 m cos X m £) X(v / 3A m ), (4) 
2v N — ; 

?7l=l 

N 

Om cos X m £ - (3 m sin \ m i) J(Vs\ m ) = 0, (5) 

m=l 

where X m = kom/N, I(x) = sine a; — cosx — (x/3)sinx with sincx = x~ x sinx and J[x) = 
xcosx — sinx. By looking at the neglected terms in this linearization procedure, we find 
that (jlj) and ([5]) accurately predict droplet equilibria provided that rjok^ <^ 1. To conform 
with this condition, together with our principal aim to examine the effects of small scale 

— 1/2 

roughness, we focus on substrate families with 770 <S 1 and 1 -C &o <C rj Q 1 . Hence, for 
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FIG. 3: Dynamics of droplet spreading over randomly varying substrates with 770 = 5 x 10 -4 and 
ko = 20, 30 and 40 using 20,000 samples for each substrate family, (a) Standard deviation of e as a 
function of time. The numerical experiments (solid line) are indistinguishable from theory (dashed 
line) for ko = 20. (b) Standard deviation of I as a function of time. The early-time asymptotics 
agree with numerical experiments up to t ~ O(30). For large times the solid lines asymptote at 
values predicted from the long-time analysis. 



a droplet with L = 0.5mm, a s = 15° and substrate topographies with amplitudes 0.5/xm 
(7/0 ~ 4 x 10~ 4 ), 7/0^0 < 1 for Iq > 77/im (beyond the region of wavelengths in Fig. [TJo, 
but Hitchcock et al. give materials with smaller amplitudes and larger wavelengths, e.g. 
for silica, 7/0 ~ 0.006 — 3.1/im, Iq ~ 65 — 110/zm). At these scales, slip is more important 
than intermolecular forces in controlling the spreading dynamics, as the study by Hocking 
on droplet spreading over flat substrates has demonstrated [18]. By the central limit 
theorem, we also see from (j3J) that e is approximately a normally distributed random variable. 
Contributions to its mean are of O (rjo) and originate from (a m sin X m £ + j3 m cos \ m £). If all 
equilibria are taken into account however, this quantity vanishes and contributions to the 
mean are of 0(j]q). Therefore, to fully assess the effects of spatial heterogeneities on wetting, 
the equilibria attained from the droplet dynamics need to be considered instead. 

For arbitrary t, the IDEs obtained for deterministic substrates in [3] are appropriately 
modified to model random spatial heterogeneities described by (j2J) and are linearized for 
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e <C 1. This calculation is rather involved and lengthy. The final equations are of the form, 



e + A (t) e = —^= y2 (a m sin k m £ + (3 m cos k m £) B (t, k m ) 
- ' N — ' 

} J (a m cos k m £ - (3 m sin k m £) C (t, k m ) 



N 



m=l 
N 



N , 

m=l 

where A(t), B(t, k m ) and C(t, k m ) are complicated functions of their arguments that we omit 
here for the sake of brevity. Their time-dependence enters through xq and its time derivative 
with xo satisfying, 3x'o ln[2xo/(Ae 2 )] = 27xq 6 — 1. The linearity of the equation for e also 
implies that this quantity is a normal variable for all times, whose variance may be computed 
explicitly. Figure [3^ shows plots of the standard deviation of function of time for 

rjo = 5 x 10~ 4 and ko = 20, 30 and 40. When ko = 20, the theoretically predicted curve is 
indistinguishable from the one obtained from numerical experiments, but as expected the 
agreement tends to degrade as the condition 770 k$ <C 1 no longer holds. Determining the 
time-evolution of the standard deviation of £, ag, explicitly is a formidable task due to the 
highly non-linear nature of the equation for £, but the early-time behavior can be found by 
linearizing about £ = 0. In Fig. [3b we show the evolution of <Je as computed from numerical 
experiments, together with the early-time behavior predicted by the linear theory for the 
same parameters as in Fig. [5b- There is excellent agreement up to t ~ (9(30); after that 
time the theoretically predicted variance goes to infinity since the linearized equation for £ 
predicts a Cauchy variable in the long-time limit. Comparing the time-scales over which a £ 
and at saturate, reveals that the droplet 'footprint', 2 (e + xq), approaches equilibrium over 
a shorter timescale compared to the time for i to reach equilibrium, which suggests that 
the droplet slides without spreading along the various substrate features to find the final 
equilibrium position. 

In the long-time limit, we cannot solve for £ explicitly since (J5|) is nonlinear and moreover 
it admits infinitely many solutions. However, the evolution towards equilibrium fixes the 
solution to (j3J) to be the stable equilibrium that is closest to £ = 0, a problem which 
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is reminiscent of the highly nontrivial 'first-passage problem' in probability theory 
Interestingly, Fig. reveals that the probability density of £ as t — >• 00, pg, is far from being 
a normally distributed random variable as the comparison with the density of the normal 
variable with the same variance demonstrates. By taking into account the mean distance 
between zeros of (jSJ) [3] , together with the fact that on average half of the closest equilibria 
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FIG. 4: Statistics of £ for ko ranging from 10 to 40 with ijo = 5x 10 _4 ,(+) and 10 -3 , (x), taken 
over 20,000 samples and N = 1,000. (a) Probability density function of £ when rjo = 5 x 10 -4 
and ko = 20 compared with a normal density of the same variance (dashed line); (b) Standard 
deviation of I as a function of ko, illustrating the excellent agreement of ([6]) with numerics. 

are unstable, we find that 



<j- f 



yk-o 2 1 1 



| sine 2^0) +C(fco 4 ) 



(6) 



in the limit ko ^> 1. It is clear that as the substrate features become rougher the droplet 
has a tendency to shift/slide less along the substrate. This behavior is confirmed in Fig. H)d, 
where we plot the theoretically predicted <jg as a function of k together with numerical 
experiments for different substrate families, confirming also the independence of I on 770- 

From our numerical experiments we also found that (e) < in the long-time limit, thus 
suggesting that surface roughness reduces wetting. Such behavior appears to contradict 
Wenzel's theory but it signifies the fact that the droplet has to overcome the energy bar- 
riers that separate the multiple equilibrium droplet states. This effect is demonstrated in 
the recent experiments of Chung et al [lj]], where spreading perpendicular to the grooves 
of parallel-grooved substrates violates Wenzel's law, and further supported by the work 
of Cox s| on wedge equilibria over (deterministic) three-dimensional rough substrates, who 
postulated that roughness- induced wetting enhancement is due to a higher-order effect which 
manifests itself when spreading does not occur perpendicular to the substrate grooves. A 
semi-analytical expression for (e) can be obtained by noting that from our numerical exper- 
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FIG. 5: Statistics of e for the same parameters as in Fig. |4) (a) Comparison of the numerically 
determined (e) and (e) approx as a function of ko; (b) Comparison of the theoretical and numerical 
standard deviation of e as a function of ko- The agreement is excellent for r/o^o ^ 1- 

iments, (a m sin \ m i + (3 m cos \ m £) = F (ko) l 7(v / 3A m )/(/c V / 3), where F appears to depend 
weakly on rjo and ko and equals to 3 for rjok^ <C 1. Based on this observation, (e) is found 
to be 

(^approx « (2 - COS 2^) + O {VOK 1 ) , (7) 

suggesting that (e) has an oscillatory behavior as fco varies. Furthermore, the mean apparent 
contact angle is now larger than the static contact angle by an amount 2| (e) \/y/3 to leading 
order in e. Figure [5h depicts a plot of (GO) as a function of ko together with the mean obtained 
in numerical simulations for substrates with 770 = 10~ 3 , considering 20,000 samples from each 
substrate family. For smaller ko, the agreement between the semi-analytic approximation 
and the numerical experiments is evident, but as the substrate becomes more rough so 
that 770/co C 1 is no longer valid and nonlinear effects become appreciable, there is a clear 
deviation towards a progressive reduction of the mean droplet radius. This implies that the 
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apparent contact angle tends to increase with substrate roughness, thus pointing towards a 
substrate-induced, hysteresis-like effect. 

The variance of e can be deduced from (j3J) by converting the Riemann sum into an integral 
by taking N —* oo, 



when k ^> 1. The theoretically predicted a e is in very good agreement with the simulated 
one as shown in Fig. [5b, where we plot a £ as a function of k when r/ = 5 x 10" 4 and 
i]q = 10~ 3 . Had we used a self-affine substrate representation, the results are qualitatively 
the same since the leading order variance of e differs from the leading order of (|SJ) only 
by a factor 3(2 - D)(g 4 - 2D - g 2 )/\(D - 1)(1 - g 4 ' 20 )], that depends on two additional 
parameters, namely the fractal dimension D and on the lower to upper cutoff wavenumber 
ratio, g. Numerical studies of o~t also confirm qualitative agreement. 

To conclude, we have presented the first detailed and systematic investigation of the 
motion of 2D droplet fronts over randomly varying shallow substrates by using a model 
derived first principles, i.e. from the governing hydrodynamic equations in the limit of 
small contact angles. Analytical predictions made for the equilibria of the two quantities of 
interest, the contact line fluctuation e and droplet shift along the substrate I exhibit excellent 
agreement with numerical experiments. However, our analysis revealed that considering 
droplet equilibria alone cannot fully determine the effects of random substrates on wetting; 
instead, their stability from the dynamical spreading problem needs to be taken into account 
as well. For arbitrary times, examination of the evolution of I and e suggests that on average 
the droplet has the tendency to slide without spreading along the substrate before reaching 
equilibrium. In the long-time limit, £ and e scale with a 2 ~ r] 2 k 2 /8 and erf ~ 57T 2 A^~ 2 /6, 
respectively. We believe that these results will motivate further analytical and experimental 
studies on the role of heterogeneities on wetting hydrodynamics. 

We thank the anonymous referee for suggesting the statistical self-affine substrate 
representation. We acknowledge financial support from EPSRC Platform Grant No. 
EP/E046029. 
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